ANL-HEP-CP-94-53 
ITP-SB-94-40 



PRINCIPAL- VALUE RESUMMATION 
FOR DILEPTON PRODUCTION 

■ HARRY CONTOPANAGOS 

■ High Energy Physics Division, Argonne National Laboratory, 

Argonne, IL 60439-4815, USA 

wo: 

^ . and 

<: 

LYNDON ALVERO, GEORGE STERMAN 

\ Institute for Theoretical Physics, SUNY Stony Brook, 

Stony Brook, NY 11794-3840, USA 

> : 

m . 

C\ ■ ABSTRACT 

m ■ 
oo ' 

We present a new method of resummation of QCD threshold distributions 

to the hard scattering function of the dilepton production cross section in 

hadron-hadron collisions. Our formulation by-passes the infrared singularities 

of the QCD running coupling through a principal value prescription, and does 

not require an explicit infrared cut-off. The resulting large corrections expa- 
ts 

nentiate, and we discuss asymptotic properties of the exponent. We present 
' a closed analytical formula for the resummed hard part in momentum space, 

and give predictions, at fixed-target energies, for the resulting cross section 
^ . which is a bounded function of its kinematic variables in the entire partonic 

phase space. 

1. Introduction 

It is well-known that perturbative QCD (pQCD) corrections to certain hadron- 
hadron scattering cross sections can be numerically important even at quite high en- 
ergies. One source of large corrections is soft-gluon radiation near partonic threshold. 
The best-studied example of this mechanism is the dilepton-production cross section, 
which is not only of great phenomenological importance but is also a true laboratory for 
testing various ideas on pQCD. As we shall see, at fixed-target energies, the resummed 
series can differ significantly from even two-loop expansions. Indeed, the perturbative 
series evaluated order-by-order does not appear to converge, and in fact, due to the 
presence of infrared (IR) renormalons, is at best only asymptotic. This suggests that 
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predictions of the normalization of such cross sections on the basis of low-order calcu- 
lations combined with parton distributions determined from deeply inelastic scattering 
are unreliable. 

Resummation formulas typically organize the large perturbative corrections in an 
exponential containing integrals of the QCD running coupling a s (fi 2 ) over various mo- 
mentum scales. The resulting expressions, however, suffer from divergences 1 associated 
with the singularity in the perturbative coupling at Aq CD . One way to deal with this 
problem is to introduce an explicit infrared cutoff. The dependence on such a cutoff is 
higher-twist 1 . The resulting dependence on the cutoff, however, is difficult to analyze, 
and is not negligible in general. 

Here we summarize recent work 2,3 ' 4,5 that advances a new resummation method 
of large threshold corrections, which we have called principal-value resummation. In 
this approach, the integral over the momentum scales of the running coupling in the 
resummed exponent is defined as a principal value integral, thus by-passing the IR 
singularities of the coupling. It can therefore be taken as a mathematically unambigu- 
ous definition of the resummed perturbative series. This then allows us to define and 
explore the structure of its higher-twist corrections, since they are defined relative to 
a well-defined perturbative expression. Furthermore, the new resummed exponent ex- 
hibits the behavior anticipated in pQCD: at each point in partonic phase space it is 
best-approximated by a particular polynomial in a s . It also becomes a purely non- 
perturbative, higher-twist expression in the non-perturbative region of the partonic 
phase space. The resulting hard part is calculable numerically as well as analytically 
(including inversion of the Mellin transform). We will summarize all of the above prop- 
erties within the context of the Drell-Yan process. 

2. Principal- Value Resummation 

We will focus on the reaction 

h 1 ( Pl ) + h 2 (p 2 )^ll(Q)+X, (1) 

where in the following we use the notation s = (pi +P2) 2 , and z = T/x a Xb, with x a , x\, 
the momentum fractions of the scattered partons. We will concentrate on the diagonal- 
flavor quark-antiquark subprocess which is singular in the soft-gluon limit z — > 1 at 
fixed order of pQCD. 

The corresponding cross section may be written as 

^ = <(g 2 )^(r,g 2 ), (2) 

where is the point-like tree-level cross section for qq — > V* — > 11, containing all 
electroweak parameters and physical dimensions and 

W v (r,Q 2 ) = j:Cj / 1 ^^F //fcl (x ,Q 2 )F //fca (x 6 ,Q 2 )a; // (r/x oa ; 6 ,a), (3) 

f -'0 X a Xfj 

where a = a s (Q 2 )/ir and Cj , in units of |e|, is the coupling constant of parton / to 
vector boson V. Specifically, Cj = e 2 and Cf — 1 + (1 — 4|e/| sin 2 ^) 2 . The sum in 



eq. ([|) runs over all relevant quark and antiquark flavors and the parton distributions 
Ff/ hl , Ff/ h2 must be in DIS scheme form. 

For the Mellin transform of the hard part, O0ff(z, a), we obtain the principal-value 
resummation formula in moment space 3 

u ff {n,a) = A{a)e E{n ' a \ (4) 

where the function A(a) exponentiates large "constants" 6 , i.e., terms that depend only 
on a and not on the parton momenta, and the resummation exponent E(n, a) depends 
on the soft gluon momenta through the correspondence n <-> 1/(1— z). It is this quantity 
that we will define in principal-value resummation. 
We may write 1 
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E(n, Q 2 ) = -l_ dC \ 7 r /7 (l - C, Q 2 ), (6) 
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and define the principal-value exponent as 3 
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where P = | x C + U | x C_, with C + (C_) any contour in the ( plane connecting 
and 1 above (below) the real axis and with a(XQ 2 ) truncated to only those terms that 
create the large threshold corrections, 

f\rfl\- f\\ a a 2 b 3 \n(l + ab 2 \n\) 

° {XQ } = a(A) = l + ab 2 lnX ~ -bT il + abtlnX)" (?) 

where b 2 = 7^(33 — 2nf) and 63 = ^(306 — 38n/), n/ being the number of active 
flavors. The functions g±, i — 1, 2 in eq. (|5]) are simple expansions of the form g^(a(X)) = 

J2"jL\ °^{^)9i an d the corresponding coefficients are factorization scheme- and process- 
dependent and can be fitted to finite-order results. What makes the procedure finite and 
useful is that we only need three of these coefficients to account for all large threshold 
corrections to all orders. These, in the DIS scheme, are 



g? = 2C F , = -2c>, g? = C t 
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It can be shown 4 that the purely perturbative regime is defined, in moment space, by 
the relation 

1 < Inn < — '— <->• 1 < n < %, (9) 
2ab 2 A 

where, in the second relation, we have used for simplicity a 1-loop relation between a 
and A = Aqcd- I n this perturbative range, a small-a asymptotic approximation leads 
to the perturbative formula for the exponent 

2 N,+l p+1 7 Ni+1 p 

E{n,a)~E(n,a,Af) = Yl E ^T. S U^^ + H E ^E^W^n, (10) 

i=l p=l j=0 i=3 p=l j=0 



where the {sj tP (i)} are calculable numerical coefficients and M = M{Q) = {iVj}, the 
set of optimum numbers of perturbative terms, is determined by optimization of the 
perturbative asymptotic approximation, eqs. @, fllPl). At typical fixed-target dilepton 
mass scales, say Q = 10 GeV, we find 3 ' 4 Af = {2, 6, 6, 2, 6, 1, 6}. 

It can similarly be shown 4 that the purely non-perturbative regime is defined 
through the relation 

1 Q 2 

Inn > — — <-» n > — . (11) 
ab 2 A 2 

Exact analytical results for the exponent in this region have been obtained 4 through a 
large-a asymptotic approximation, but suffice it to say that for extremely large values 
of n, 

(l) (l) (Ik 

E(n,a) ~ -^-ln21nn- ^-\n(ab 2 \nn) - ^-^\n 2 (ab 2 \nn). (12) 
v ' b 2 b 2 Abl 

This shows that the exponent tends to — oo at the edge of phase space z ~ 1 {i.e., 
n — > oo). As a result, the corresponding resummed cross section is finite, because the 
hard part is exponentially damped at that threshold. This results in the principal-value 
cross section being a bounded function of r throughout its kinematic range, in contrast 
to the finite-order cross section that diverges as a power of the logarithm as r — ► 1. 

To obtain the resummed hard part in momentum space, we perform an inverse 
Mellin tranform on eq. (^) analytically 2 . The final result for the resummed hadronic 
function, eq. (^), is 

W v (r,Q 2 ) = J2CjA(a) £ dz[l + H(z, (^7^) > (13) 

where the parton flux is defined as 

Ffjirlz) = !\{l-r/{zx a x b ))^^F f , hl {x a ,Q 2 )F f/h% {x b ,Q 2 ), (14) 

JO X a Xjy 

and the hard function TC is given by 2 ' 4 

a) s [ iz 'w^7) T i 1 + Pi (rb- a )) sm h (rb- a )) • (15) 

where T is the Euler Gamma function and P\(x, a) = -g^-E(x, a). This hard function 
may be used for predicting both the mass-distribution for the process and its derivative 
with respect to rapidity or xp 5 . 

3. Comparison with experiment and conlusions 

Because of space limitation, we will treat only a single example: the E537 fixed 
target experiment at Fermilab 7 , which observed Drell-Yan pairs produced by antipro- 
tons on a tungsten target at y/s = 15.3 GeV (V — 7 in eqs. @ and @). In fig. 1, 
we plot the data as well as the one-, two-loop and resummed predictions for the mass 



distribution Q 3 dcr/dQ integrated over xp > 0. We used a parton distribution 8 based on 
a fit to the purely DIS data of CDHS. Notice that agreement between the resummation 
prediction and the data is acceptable, and better than that of the finite-order predic- 
tion. Parton distributions with larger sea-quark content tend to give larger predictions 
in the resummed cross section. 

We conclude that principal-value resummation is an appropriate way to organize 
consistently large threshold corrections in pQCD to all orders in perturbation theory, 
and that experimental results can be successfully described in this framework. This 
also suggests that whatever higher-twist part of the theory is not contained within 
pQCD may be supressed, even in the domain of phase space probed by fixed-target 
experiments at momenta like those cited in fig. 1. The inclusion of the information 
contained in resummation formulas may also be important for the determination of 
parton distributions by the combination of DIS and hadron-hadron data. 




Fig. 1. E537 experiment. Short-dash: 1-loop; Long-dash: 2-loop; Solid: Resummed. 
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